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It is known that the anomalous Hall conductivity (AHC) in a disordered two dimensional electron 
system with Rashba spin-orbit interaction and finite ferromagnetic spin-exchange energy is zero in 
the metallic weak-scattering regime because of the exact cancellation of the bare-bubble contribution 
by the vertex correction. We study the effect of inhomogeneous longitudinal electric field on the 
AHC in such a system. We predict that AHC increases from zero (at zero wavenumber), forms a 
peak, and then decreases as the wavenumber for the variation of electric field increases. The peak- 
value of AHC is as high as the bare-buble contribution. We find that the wave number, q, at which 
the peaks occur is the inverse of the geometric mean of the mean free path of an electron and the 
spin-exchange length scale. Although the Rashba energy is responsible for the peak- value of AHC, 
the peak position is independent of it. 
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I. INTRODUCTION 



Hall resistance in a ferromagnetic material is known to 
have two contributions: 1) proportional to magnetic field, 
2) proportional to magnetization. While the former is re- 
sponsible for conventional Hall effect, the latter is due to 
anomalous Hall effect ( AHE) . The anomalous Hall coeffi- 
cient is the zero field Hall resistance which is estimated by 
extrapolating Hall resistance up to zero magnetic field. 
The physical mechanisms which contribute to the anoma- 
lous Hall conductivity (AHC), are broadly categorized 
into two classes: intrinsic and extrinsic. While the prop- 
erties of the bands are responsible for the former, the spin 
dependent scattering of electrons belong to the extrinsic 
mechanism. Karplaus and Luttingcr 1 proposed intrinsic 
mechanism due to the spin-orbit interaction (SOI) for 
the polarized electrons in ferromagnetic materials. This 
mechanism was later realized as the outcome of the Berry 
phase 2 from band topology. Smit 3 proposed an extrinsic 
mechanism which is known as 'skew-scattering' (SS) in 
which electrons scatter from impurities with spin asym- 
metry. In addition, contribution due to side-jump (SJ) 
scattering was put forward by Bcrgcr 4 through a semi- 
classical theory. 

The debate on the true mechanism for AHE grew 
stronger with the contradicting results 5 for various sys- 
tems. The linear transport theories, which are now 
mostly used to study AHE, seem to have settled some 
of the confusing issues in the recent years. Although, 
different approaches of this theory often tend to conflict 
with each other, there are some successful attempts 6,7 to 
establish links between them. To this end, the two di- 
mensional ferromagnetic- Rashba gas, i.e., electrons with 
Rashba SOI and ferromagnetic exchange coupling has 
proved to be a simple and elegant model to work with 
because it captures almost all the important and fun- 
damental features relevant to AHE. Due to the Rashba 



SOI, characterized by the Rashba energy A R , degeneracy 
is lifted by the formation of two chiral subbands with 
opposite helicity. Further, ferromagnetic spin-exchange 
energy produces a gap between these subbands. 

A number of studies on AHC have been performed 
within the model of two-dimensional ferromagnetic- 
Rashba gas, using the Kubo formalism 6 ' 8 ~ 14 ,the Keldysh 
technique 6 ' 15 > 16 , and the semiclassical treatment of 
Boltzmann transport equation 17 ' 18 . According to the 
studies in Kubo formalism, the magnitude of AHC is 
sensitive to whether only one or both the subbands are 
occupied by the electrons, i.e., whether or not the Fermi 
level is situated below the band-gap created by exchange 
energy, A ex . In the presence of non-magnetic impurities, 
both the Fermi sea and Fermi surface contributions in 
AHC vanish 11 when e F > A ex , where e F is the Fermi 
energy. The bare-buble contribution is exactly canceled 
by the vertex correction 10 ' 11 due to SJ mechanism, and 
the SS term is absent 12 . On the contrary, all these terms 
are finite for e F < A ex . A comprehensive analysis of 
these works can be found in the detailed study of Nun- 
ner et al 11 . Using the Keldysh technique in a T- matrix 
approximation, Onoda et al 15 ' 19 , however, find a strong 
SS contribution that leads to non-zero AHC even when 
e F > A ex . They also observe that the AHC changes 
sign when the Fermi level enters into the band gap. The 
Keldysh technique automatically incorporates the higher 
order Born scatterings and enables them to arrive at such 
striking conclusions when the Fermi energy is located 
around the anticrossing of band dispersions in momen- 
tum space. Kovalev et al 16 later included the higher order 
scatterings explicitly in the Kubo formalism to match the 
results of Onoda et al 15 and were able to unify the these 
approaches at least for a certain range of parameters. 

There is no controversy, however, in the metallic weak- 
scattering regime (e F r > 1), where r is the mean scat- 
tering time of an electron. When the Fermi level goes 
far above the band gap (e F > A ex ), i.e., when both the 
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chiral subbands are occupied, the AHC vanishes 5 ' 7,11,20 , 
making this regime a somewhat uninteresting one. We 
however note that AHE is intimately related to spin hall 
effect 21,22 and there are instances of enhancement of spin- 
Hall conductivity 23 and occurrence of a induced trans- 
verse force 24 due to inhomogeneity of the applied electric 
field. In this paper, we study the effect of inhomogeneous 
longitudinal electric field (E q || q) on the AHC in the 
metallic weak-scattering regime when both the subbands 
are occupied. 

We find, within the Kubo formalism, that AHC, a AH , 
is zero at q = 0, increases with q and forms a peak at cer- 
tain q before it decreases to zero at large q. a AH depends 
on the parameters e F r, A ex r, and Art. We find that the 
peak-value of AHC is proportional to (Art) 2 (c f t)^ 1 and 
the peak occurs when ql e ff ~ 1, where l e jf = (llex) 1 ^ 2 , I 
is the mean free path of an electron, and l ex = v F / (AA ex ) 
is the length-scale corresponding to spin-exchange, with 
v F being the Fermi velocity. We note that the momen- 
tum at which the peaks occur does not depend on A R . 

The paper is organized as follows. The following sec- 
tion contains a brief formalism for the evaluation of the 
anomalous Hall conductivity. In Section III, we present 
our results and discuss about the relevant scales and lim- 
its. Finally, we summarize our results in Section IV. 



II. ANOMALOUS HALL CONDUCTIVITY 

A. The model 

A system of spin-polarized two dimensional electron 
gas with Rashba spin-orbit interaction in a disordered 
environment may be expressed by the Hamiltonian 



H = I — + V(r) I a +a(ax V)-z- A ex a z , (1) 



where m is the the effective mass of an electron, a is 
the Rashba spin-orbit coupling parameter, A ex is the ex- 
change energy which favors one kind of spin over the 
other, a represents three Pauli matrices and <7o is 2 x 2 
unit matrix. Here V(r) is the spin independent disor- 
der potential for the randomly located 5-function impu- 
rities and has the form V(r) = ViS(r — Hi) satisfying 
(Vi) = , (V 2 ) =V 2 ^0 and (V?) = V 3 £ 0. (We have 
set the unit h = 1 and c = 1.) 

In the absence of disorder, the above Hamiltonian can 
be exactly solved with the eigen values 



(2) 



where k is the momentum of an electron, s = ± is the 
label for two chiral subbands produced due to Rashba 
spin-orbit interaction and Cfc = y A 2 X + a 2 k 2 . The cor- 
responding Fermi momenta K s are given by 



K 2 = 2m [e F - < s ] 



(3) 



where e F is the Fermi energy and 



C, = J A 2 ex + A 2 R + m 2 a 4 - s ma 2 , 



(4) 



with An = ay'2me F = ak F . Here k F is the mean Fermi 

k 2 

momentum, given by e F = j 2- . The density of states at 
the Fermi level for the two subbannds are 



A 2 R + m 2 a 4 



(5) 



where vq is the density of states for each spin in pure two 
dimensional electron gas. The free retarded (advanced) 
Green's function of the system can then be expressed as 



G 



ra, ,_lr^ op + sa(cr x k) ■ z - A ex cr z )/Ck 



(0) 



where = — ep. The self energy for the electrons 
due to scattering is then found to be 



^R(A) 



4r^ 



T- 



(7) 

in Born approximation, where r = (27r^onjV2) 1 is the 
mean scattering time for an electron, and rii represents 
the concentration of impurity. We thus find the disorder- 
averaged Greens functions as 



1 (T + sa(o- x k) • z - A ex a z )/( k 



2 ti 



(8) 
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Here and Ck get renor- 



malized to and £fc due to the presence of A ex in- 
stead of A ex in their respective expressions. When both 
the helicity subbands are occupied, v+ + V- = 2vq and 
= V— /C- and therefore exchange energy does not 
get renormalized and e becomes e + i/2r. 



B. Kubo Formula 

Using Kubo formula, anomalous Hall conductivity may 
be expressed as 



T AH ~- 

'yx 



1 



2tt J (2tt) : 



Tr 



{I (k+|)+f a (q)}^(0) 



+ ~1 G^ +q (0) 



(9) 



where charge current operators j x (k) = e(^-<7o — aa y ) 
and j y (k) = e(-^-cr + cut x ), and T x (q) represents the 
correction to the vertex of j x due to both side-jump and 
skew scattering contributions (see fig.l) which can be cal- 
culated by solving the self-consistent equation 
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FIG. 1. (Color online) Diagrammatic representation of the 
self-consistent equation [10] including side-jump (ladder) and 
skew-scattering corrections. The small wavy lines with a point 
vertex represent the bare charge current vertex j x , and the 
small wave lines with dressed vertex represent T x . Retarded 
and advanced Green's functions are represented by the solid 
lines with arrows in the opposite directions. Dashed lines refer 
to the impurity scattering. Lowest order side-jump (a) and 
skew-scattering (b,c) diagrams; diagrams (d,e,f) represent 
vertex correction to the diagrams (a,b,c) respectively. Here 
vertex correction due to both side-jump and skew-scattering. 



C. Numerical evaluation 

The numerical evaluation of a^ 1 via Eqs. (11) and 
(12) proceeds with the same spirit of Refs. 14 and 
23. Expanding J x (q) in the Pauli matrix basis: J x = 
J2 a =o Ja&a, Eq.(12) can be rewritten as 



xG£(0) = 



Equating the coefficients of the four Pauli matrices in the 
above equation, we find a 4 x 4 matrix equation: 



1 - moo -mu -77122 -TO33 

— moi 1 - toio im 23 -im 32 

-m 2 -imi 3 1 - to 20 im 3 i 

-m 3 im 12 -im 2 i 1 - m 30 
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(14) 



where m a p (q) and zp (q) are given by 
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f*(q) =7i 



72 



72 



dk' - 

(2^ G ^ +q (0) 
(2tt) 2 k ' +q 



Jx 

dk" 



(27T) 



2 ^k"+q(0) 



> Ik' 

aK Ail 
(2tt) 2 k " +q 

r-f^q)] G&,(0) 

q 



+ f x (q) G&(0) 



k" 



+ P,(q) 



dk' 



G£(0), 



where, 71 = riiVi and 72 = rnV 3 . It is established from 
a number of studies 6 ' n > 12 that the skew-scattering con- 
tribution to the homogeneous (q = 0) anomalous Hall 
conductivity vanishes provided both the chiral subbands 
are occupied. This, however, is true for any q since the 
last two terms in Eq.(10) cancel each other because of the 
identity J G^' A = ^mv Q a . Equation (9) therefore 
reduces to 



rr AH =— 

y 2ir 



dk 

{3* ( k - 



Tr 

q 

" 2 



f ) ^ +q (0) 



>y { k 
Jx(q)}^(0) 



(11) 



with J x (q) containing side-jump types of scattering con- 
tribution of (q) only, i.e., 

J -M = ^/<^w°>{Mt+f) 

+j,(q)}6£(0). (12) 
We numerically evaluate J x (q) and hence below. 
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G* +n a a G£a a (15) 



2/? (q) = 



4irv T 



■Tr 



^2^ +q (o)3, (k + %)6£(o)* p 



(27T) 



(16) 

Inverting the matrix equation (14), we find J Q and then 
substituting these in Eq.(ll), we obtain a^ 1 . 



(10) 



III. RESULTS AND DISCUSSIONS 



In the metallic regime (e F r 3> 1), there is no net cr^f* 
when the applied electric field is uniform. When both 
subbands arc partially occupied, the bare bubble contri- 
bution 



A 2 



+ (27) 2 



(17) 



is exactly canceled by the ladder vertex correction, a AH , 
which is numerically equal to a AH but opposite in sign. 
In the super clean limit (r — > 00), the above value agrees 
with the previous result . 

In this paper, we are considering the effect of inho- 
mogencous longitudinal electric field, i.e., q || Eq on 
the Anomalous Hall conductivity. We define three di- 
mensionless parameters out of four energy scales e F , A/?, 
A ex and 1/r. They are e F r, Art and A ex T. We numeri- 
cally calculate a^ 1 for different sets of these parameters 
given in Tablc-1, using Eqs. (11), (14)-(16) when the 
electric field is assumed to be applied along x-axis, i.e., 
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FIG. 2. (Color online) (a) The variation of anomalous Hall 
conductivity in the unit of e 2 /(2ir) with q x /k F for different 
sets of parameters (i)-(v) shown in Table- 1. The peaks in a AH 
occur at different valued of q x /k F depending on the choice of 
the parameters. Scaling q x /k F by (e F r) _1 (A ea; r) _1 ' /2 , a AH 
is replotted (b) against qml e tt f° r all the chosen sets of pa- 
rameters. Peaks in (b) occur when qxl^tf ~ 1. 



Set 




A ex T 




-a~r[e</(2n)] 


(i) 


0.2 


0.4 


fO.O 


0.0018 


(ii) 


0.2 


0.8 


5.0 


0.0034 


(hi) 


0.25 


0.4 


5.0 


0.0061 


(iv) 


0.2 


0.2 


2.5 


0.0071 


(v) 


0.4 


0.4 


5.0 


0.0112 



TABLE I. Five sets (i)-(v) of three parameters Art, A ex T, 
and € f t used to calculate a AH for the Fig. (2). The last 
column shows the value of a AH calculated using Eq.(17) in 
the unit of e 2 /(2tt) for these sets of parameters. 



As evident in Fig. 2(a), the magnitude of a in- 
creases, forms a peak, and then decreases with the in- 
crease of q x . We note that the peak values for the differ- 
ent sets of parameters are close to the corresponding bare 
contributions a AH listed in Table- 1. This suggests that 
the effect of inhomogcneous electric field on the AHC 
is substantial. The peak positions and the maximum 
values of a AH depend on these three parameters. How- 
ever, if we scale q x in the unit of inverse effective mean 
where I — v f t is the mean free 
v F /(AA ex ) is the length 



free path, l e ff = y 
path of an electron and l e 
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FIG. 3. 



(Color online) Contour plot of scaled anomalous 
AH (e F T)(A R T)~ 2 versus jp- and 



Hall conductivity, a AH 



The darkest linear patch, where the magnitude 



of a is maximum, occurs when qxl ell ~ 1. 



scale corresponding to the exchange energy, the peak po- 
sitions occur at q x leff ~ 1; irrespective of the values 
of the parameters (see Fig. 2(b)). On the other hand, 
the magnitudes of the peaks depend on two parameters: 
(A fl r) 2 (e F r) _1 . Figure 3 shows a contour plot 
2 as a function of and 

{A ex T) 1 / 2 {e F Ty 1 , i.e., (fc F ; e// ) _1 . Note that for the 
nearly-linear darkest patch, in which the magnitude of 
a AH is maximum, corresponds to q x l e ff ~ 1- 

These results suggest that if an electric field of the 
form E x = Eoe iqx is applied along the x-axis in the 
plane of a two dimensional electron gas, the AHC be- 
comes finite. As the maximum magnitude of AHC occurs 



of & AH = a AH (e F r)(A RT ) 



at q x /kp = (A ex r) 1 / 2 (e F r)~ 1 , the periodic variation of 
the applied field should be properly tuned in accordance 
with the relevant parameters of the system to increase the 
transverse charge current. Clearly, the nature of this pe- 
riodicity depends on magnetization, disorder and carrier 
concentration (i.e. on A ex ,r, and e F respectively) and 
not on Rashba parameter a. This is illustrated in Fig. 
4. An increased variation of the electric field is necessary 
to produce maximum AHC in a system with higher leff, 
i.e., with increasing magnetization and disorder, and de- 
creasing electron density. 



IV. SUMMARY 

To summarize, we have calculated the anomalous Hall 
conductivity for a two dimensional ferromagnetic Rashba 
system with an electric field inhomogcneous in nature. 
It vanishes for a homogeneous electric field in the weak 
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I I M" «; — I ■ ■ ■ ■ 

FIG. 4. (Color online) Schematic variation of the inhomo- 
geneous electric field E x , causing resonance in AHC : The 
periodicity is illustrated by the color gradient. For increased 
magnetization and disorder and decreased Fermi energy, i.e.,, 
for increased Z e //, the spatial variation of the electric field 
should be increased to obtain this resonance. 



scattering regime when both the subbands are occupied. 
We have shown that the anomalous Hall current can in- 
crease considerably due the variation of the electric field 
and this comes from the bare bubble and ladder vertex 
correction and not from the skew scattering. Dealing 
with a realistic parameter range, we have found that the 
length scale associated with the peaks in the conductivity 
is proportional to the geometric mean of the two length 
scales arising from disorder and the exchange interaction. 
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